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Pattern formation in wet granular matter under vertical vibrations 
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Experiments on a thin layer of cohesive wet granular matter under vertical vibrations reveal 
kink separated domains that collide with the container at different phases. Due to the strong 
cohesion arising from the formation of liquid bridges between adjacent particles, the domains move 
collectively upon vibrations. Depending on the periodicity of this collective motion, the kink fronts 
may propagate, couple with each other and form rotating spiral patterns in the case of period 
tripling, or stay as standing wave patterns in the case of period doubling. Moreover, both patterns 
may coexist with granular ‘gas bubbles’ - phase separation into a liquidlike and a gaslike state. 
Stability diagrams for the instabilities measured with various granular layer mass m and container 
height H are presented. The onsets for both types of patterns and their dependency on m and 
H can be quantitatively captured with a model considering the granular layer as a single particle 
colliding completely inelastically with the container. 

PACS numbers: 45.70.Qj, 45.70.-n, 45.70.Mg 


I. INTRODUCTION 


From Chladni figures p, 
dune formation in nature 
lar mixtures in industries 


to Faraday heaping Q, from 
to the segregation of granu- 
, pattern formation in gran¬ 
ular matter has been attracting interest from physicists, 
engineers and geologists over centuries IMS- In the 
past decades, many intriguing instabilities have been 
discovered, such as triggered avalanches [H, propagat- 
ing fronts 0, heap corrugation and transport [10, O , 
convection fingering in air [13 or in water [T^, 

and stratification patterns of granular mixtures driven by 
avalanches [a or by horizontal vibrations HEl- Par¬ 
ticularly, standing wave patterns in a vertically agitated 
dry granular medium have been investigated extensively, 
both in quasi-two-dimensional and three-dimensional 
systems [^, . Depending on the vibration strength, 

the system undergoes a period doubling cascade accom¬ 
panied by stripe, square, hexagonal and spiral patterns 
as well as kink waves . The dispersion relation of 

the subharmonic patterns was found to be reminiscent of 
gravity waves in a fluid [2^ , suggesting the possibility to 
describe vibrofiuidized granular matter as a continuum. 
Localized excitations, also coined as oscillons, were also 
discovered in the hysteresis region of the primary square 
or stripe pattern forming instabilities [29|. 

Due to the strong cohesion arising from the formation 
of liquid bridges between adjacent particles, partially wet 
granular matter exhibits a different pattern forming sce¬ 
nario: Three armed rotating spiral pattern, which mani¬ 
fests a peculiar period tripling bifurcation, was found to 
dominate [sO]- The spiral arms correspond to the kinks 
separating domains that collide collectively with the con¬ 
tainer at different phases. Further characterizations on 
the dynamics of spiral arms revealed that the rotation 
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frequency is finite at the threshold and grows linearly 
with the peak vibrating acceleration [3l| . 

As the granular layer moves collectively within the con¬ 
tainer, any other periodicity beside 3 may well exist in 
the system. Thus, it is intuitive to seek for patterns 
with various periodicities so as to find a clue on what 
determines the periodicity and why the period tripling 
state is preferred. Here we present the stability diagrams 
measured with a systematic variation of two associated 
parameters: Granular layer mass m and container height 
H. We find, in addition to period tripling patterns, pe¬ 
riod doubling standing wave patterns as well as phase 
separations [^, . Both period doubling and tripling 

regions found in the experiments match the predictions 
of a single particle bouncing model considering collisions 
with both the bottom and the lid of the container, illus¬ 
trating the roles that the two parameters play in deter¬ 
mining the periodicity. Thus, this investigation provides 
possible pathways to control the periodicity of patterns 
in vertically agitated wet granular matter. 


II. EXPERIMENTAL SET-UP 


Figure [T] shows a schematic of the experimental set¬ 
up. Cleaned glass beads (SiLiBeads S) with a di¬ 
ameter of d = 0.78 mm and 10 % polydispersity, af¬ 
ter being mixed with purified water (specific resistance 
18.2MU-cm, LaborStar TWF), are added into a cylin¬ 
drical container with a fixed mass of 768 g, an inner ra¬ 
dius R = 8 cm and a height H. The filling fraction 
of the particles is defined as ^ = m/p^H), where 
Pg = 2.50 g • cm“^ is the density of the glass beads. 
The height of the granular sample is estimated with h ~ 
m/{nR^p^T]), where p corresponds to the packing den¬ 
sity. Here the density of random close packing p = 0.64 
is used as an approximation. The sample is kept within 
the pendular state (i.e., cohesion arising mainly from liq¬ 
uid bridges formed between adjacent particles at contact) 
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FIG. 1. (Color online) Schematic of the experimental set¬ 
up with definitions of the granular layer thickness h and the 
container height H. The magnihed view of the wet granular 
sample illustrates a liquid bridge formed between two adjacent 
particles, zo, /, and t denote maximum vibration amplitude, 
frequency, and time in the laboratory frame correspondingly. 


through keeping the liquid content W = V^/Vg ^ 1.6 %, 
where and Vg are the volume of the wetting liquid and 
that of the glass beads correspondingly. The container 
is agitated vertically against gravity with an electromag¬ 
netic shaker (Tira TV50350). The frequency / and am¬ 
plitude Zo of the sinusoidal vibrations are controlled with 
a function generator (Agilent FG33220) and the dimen¬ 
sionless acceleration F = ^.'K^pzo/g is measured with an 
accelerometer (Dytran 3035B2), where g is the gravita¬ 
tional acceleration. The collective behavior of the sample 
is captured with a high speed camera (IDT MotionScope 
M3) mounted above the container. The camera is trig¬ 
gered by a synchronized multi-pulse generator to capture 
images at fixed phases of each vibration cycle. More de¬ 
tails on the set-up control can be found in Ref. j^ . 


III. STABILITY DIAGRAM 

Figure [2] shows the stability diagram measured with 
decreasing F at various /. Increasing F yields quan¬ 
titatively the same threshold for the onset of patterns. 
The boundaries between various states are determined 
through an average of three runs of experiments and the 
uncertainty corresponds to the step of F. There exists 
slight hysteresis (^ 10 %) for the critical acceleration 
at which the transition from a solidlike (i.e., no parti¬ 
cles moving) to a liquidlike (i.e., particles moving around 
while keeping contacts with their neighbors) state, as well 
as the phase separation into a coexistence between a liq¬ 
uidlike and a gaslike state (‘gas bubbles’), arises. Such a 
hysteretic behavior is in agreement with former investi¬ 
gations [32|,[^. In the gaslike state, most of the particles 
move individually without contacts to the others, and the 
granular layer expands to fill the whole container. The 
melting transition is measured through monitoring the 
changes of subsequent top view images, as the transi¬ 
tion tends to start from the free surface of the granular 
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FIG. 2. (Golor online) Stability diagram showing the col¬ 
lective behavior of the granular sample in the experiments, 
(a), (b) and (c) are typically snapshots for period doubling 
standing wave pattern, period tripling rotating spirals, and 
granular ‘gas bubbles’ (i.e., the coexistence of a liquid- and 
a gaslike state). The rotating spiral pattern shown in (b) is 
an average of three images taken at consecutive vibration cy¬ 
cles to enhance the contrast. Other parameters: m = 113 g, 
h ~ 3.5 mm, H = 10.6 mm, and ^ = 0.21. 


layer IdHj. 

If the vibration strength increases further, a standing 
wave pattern emerges at a relatively low vibration fre¬ 
quency / = 50 Hz. As the typical snapshot (a) illustrates, 
the pattern is composed of domains split by curved inter¬ 
faces that appear bright. A measurement of the surface 
profile with the laser profilometry method [36| reveals 
that the bright interfaces correspond to the kinks sepa¬ 
rating regions with different heights. At even lower fre¬ 
quencies, visual inspection is hindered by particles stick¬ 
ing on the lid of the container. In contrast to the non¬ 
propagating kinks (i.e., standing waves), a rotating spiral 
pattern is found to dominate for / > 60Hz. As shown in 
(b), the typical pattern is composed of three armed ro¬ 
tating spirals or propagating fronts. The spirals typically 
have meandering cores and may couple with each other 
through sharing arms. Different from the standing wave 
pattern described above, the kink waves propagate as pe- 
riod tripling breaks the spatiotemporal symmetry Q . 

As F increases further above the pattern forming re¬ 
gion, the system may evolve into a featureless liquidlike 
state or phase separate before entering the homogeneous 
gaslike state, depending on the vibration frequency. As 
shown in (c), the density in the gaslike region is much 
lower than the surrounding liquidlike region. Temporal 



























































































3 


1 2 3 4 48 49 



FIG. 3. Upper panel: Snapshots of a typical period doubling pattern taken at a fixed phase of consecutive vibration cycles 
(indicated as numbers) with / = 50 Hz and F ^ 17.3 g. Lower panel: Spatially resolved covariance of subsequent images 
captured as the granular layer collides with the container in another run of experiment with a closer view on the center region 
of the container. The brightness is related to the averaged mobility of the particles. 


fluctuations may occur in the liquidlike region, while no 
patterns are observed in the ‘gas bubble’ region. The 
transition from the liquidlike to the gaslike state is de¬ 
termined from a sudden increase of the noise level of the 
signal from the accelerometer, because individual colli¬ 
sions between the particles and the container are much 
more prominent in the gaslike state in comparison to the 
liquidlike one. Such an enhanced noise level also leads to 
larger error for the onset of the gaslike state in compari¬ 
son to the other state boundaries. 

Figure [3] shows the time evolution of a standing wave 
pattern. It is a labyrinth type pattern [see also Fig. UKa)] 
and different runs of experiments yield dramatically dif¬ 
ferent morphologies of the pattern. Here, it behaves like 
two fingers slightly wrapped around each other. The im¬ 
ages in the upper panel taken at a fixed phase of sub¬ 
sequent vibration cycles indicate that the pattern has a 
periodicity that doubles that of the vibrating plate. This 
is similar to the period doubling patterns observed in 
agitated dry granular matter [^. However, there are 
no periodic structures (e.g., squares or stripes) observed 
here. This difference might be attributed to the strong 
cohesion from the presence of liquid bridges: The injected 
momentum in the vibrating direction is more difficult to 
be transferred to the horizontal direction in comparison 
to the dry case. Moreover, the standing wave pattern is 
also reminiscent to the ‘phase bubbles’ observed in verti¬ 
cally vibrated dry granular layers as they are both 
composed of kink separated domains vibrating with a 
phase difference. In comparison to the transient ‘phase 
bubbles’ that shrink and disappear within tens to hun¬ 
dreds of vibration cycles after nucleation, the standing 
wave pattern observed here is more stable. Within a 
time scale of a few hundred vibration cycles, the mor¬ 
phology of the standing wave pattern stays the same. As 
time evolves up to 10^ vibration cycles, a target pattern 
sometimes arises. 

A combination of the snapshots taken after the gran¬ 
ular layer collides with the container (upper panel of 
Fig. [3|) and spatially resolved covariance C{x,y) (lower 


panel) reveals the dynamics of the pattern. C{x^y) is 
calculated with 

d 

C{x,y)= ^ /t(a:;+Aa:,y+A 2 /)/t+At(a:;+Aa:,y+Ay), 

Aa:, Ay=0 

( 1 ) 

where It{x^y) corresponds to the intensity of an image 
captured at time t. At is the time step between consecu¬ 
tive frames, and d denotes the particle diameter in pixels. 
In order to resolve the mobility of the particles, the im¬ 
ages are captured with 10 frames per vibration cycle [i.e., 
Ai = 1/(10/)]. 

In the first vibration cycle, one side of the front 
(marked with I) collides with the container and collects 
kinetic energy there. Consequently, the mobility of parti¬ 
cles in this region enhances so that an expansion into the 
neighboring region II is expected. In the next vibration 
cycle, the front will be pushed back as region II collides 
with the container. Therefore, the fronts swing back and 
forth at the time scale of vibration cycles. Such a process 
is better illustrated with (7(x, y) shown in the lower panel 
of Fig.[3l as the black (more mobile) and white (less mo¬ 
bile) regions switch at each vibration cycle. For the case 
of non-cohesive dry granular matter, the friction between 
the kink separated granular layers moving out of phase 
can induce convection in the vicinity of the kink region 
for both two- [2^ and three-dimensional [38[ systems. 
For the case of wet granular matter with strong cohesion 
between adjacent particles, it is still unclear whether con¬ 
vection can be triggered or not. This will be a topic for 
further investigations. 

As the granular layer collides with the container collec¬ 
tively, we speculate that the mass m of the granular layer 
(or corresponding layer thickness h) and the height of the 
container H play essential roles in determining the period 
for the granular layer to travel within the container. In 
order to explore such an influence, we compare stability 
diagrams measured with a systematic variation of both 
parameters. 
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FIG. 4. (Color online) Stability diagrams for various com¬ 
binations of sample mass m and container height H but the 
same filling fraction ^ = 0.18: (a) m = 96g, h ^ 3.0 mm and 
H = 10.6 mm; (b) m = 111 g, h ~ 3.5 mm and H = 12.3 mm; 
(c) m = 81.6 g, h ^ 2.5 mm and H = 9.0 mm. The bound¬ 
aries between various regimes are obtained the same way as 
in Fig. O Overlapping between various regimes corresponds 
to coexistence. 


Figure [U(a) shows the stability diagram for the same 
container height H = 10.6 mm but smaller m in compar¬ 
ison to Fig. [21 Such a change of the sample mass leads to 
the following features in the stability diagram: (i) The 
period tripling pattern region expands to / = 50 Hz in a 
r region smaller than that of the standing wave pattern. 
This suggests that the periodicity of the pattern is not 
solely determined by the vibration frequency, (ii) The 
phase separation regime expands to both lower frequency 
(/ = 60 Hz) and lower T, particularly in the high fre¬ 
quency regime. The enhanced probability for ‘gas bub¬ 


bles’ to nucleate can be attributed to the granular layer 
thickness: The thinner the granular layer, the lower the 
energy barrier for a ‘gas bubble’ to nucleate because the 
number of liquid bridges to break along the vibration di¬ 
rection is reduced, (iii) The period tripling pattern and 
‘gas bubble’ regimes can overlap, leading to a state in 
which propagating kink fronts or rotating spirals are ob¬ 
served in the liquidlike region surrounding the ‘gas bub¬ 
bles’. 

Moreover, we discuss the influence of H through a 
comparison between Fig. |4](b) and Fig. [21 as they share 
roughly the same m. As increases, the phase separa¬ 
tion region expands and the thresholds for both patterns 
grow faster with /. The latter feature suggests that the 
influence of H might be associated with the departure ve¬ 
locity r/ f) required for the granular layer to reach the 
container lid. The strong influence of H indicates that a 
collision between the granular layer and the container lid 
is essential for the observed patterns. Note that the max¬ 
imum elevation height of a particle bouncing completely 
inelastically with the container (see Sec. [TV| for details) 
is typically much larger than the container height in the 
pattern forming region. 

The stability diagrams presented in Fig. |4] have the 
same filling fraction 4> but various m. A comparison of 
(b) or (c) with (a) reveals the ‘gas bubble’ region shrinks 
(expands) as granular layer is thicker (thinner). This is 
similar to the above comparison between various m but 
the same H. For a granular layer thickness of roughly 3- 
4 particle diameters (c), the ‘gas bubble’ region expands 
further into the period doubling region. For / > 50 Hz, 
all pattern forming instabilities coexist with phase sep¬ 
aration; that is, they appear in the liquidlike phase sur¬ 
rounding the ‘gas bubbles’. Moreover, the comparison 
also reveals that the onset for both patterns grows faster 
(slower) with / as m increases (decreases), reminiscent of 
the above comparison between various H but the same 
m. The similarities to the variation of m or H individu¬ 
ally suggest that the influence of m and H is not coupled 
with each other through the filling fraction. 


IV. MODEL 

As shown in Fig. [5l the model considers the granular 
layer as a single wet particle colliding completely inelas¬ 
tically with the container. The weight of the particle is 
G = —mgez with Cz a unit vector pointing in the ^ direc¬ 
tion. The magnitude of the cohesive force is estimated 
with Fb = nNcCrd cos 0^ where Nc = is the 

number of contacts between the granular layer and the 
container, a = 0.072 N/m is the surface tension of water, 
and ^ = 0 is the contact angle. Here the area fraction 
V 2 B = 0.84 is estimated with the random close packing 
of spheres in two dimensions The total force act¬ 

ing on the particle gives rise to two critical accelerations: 
a\ = {n — l)g for the container lid and Ub = — (n + 1)^ 
for the container bottom, where n = Fb/(m 5 f) is the ra- 
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FIG. 5. (Color online) A sketch of the model considering 
the granular layer as a single particle undergoing completely 
inelastic collisions with the lid and bottom of the container. 
The red (gray) trajectories correspond to the free flying period 
of the particle. G and Fh denote the weight of the granular 
layer and the cohesive force between the granular layer and 
the container. 
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FIG. 6. (Color online) Periodicity predicted by the numerical 
analysis based on the single particle bouncing model. The 
parameters are chosen to match those in Fig. [J] 


tio between the cohesive force and the gravitational force 
acting on the particle. 

In the numerical analysis, the particle is placed initially 
on the container bottom. As the criterion for the accel¬ 
eration of the container a = —sin27r/t < Ub is 
fulfilled, the particle starts to detach from the container 
bottom. If the detachment eventually leads to a break of 
the capillary bridge, the rupture process succeeds and the 
particle starts a parabolic flight. The trajectory of the 
particle during the rupture process is calculated consider¬ 
ing the decay of the capillary force with increasing sepa¬ 
ration distance 0. From the trajectory of the parabolic 
flight and that of the container bottom or the lid, we de¬ 
termine the next collision time numerically with Brent’s 
algorithm . Considering the thickness of the granular 
layer, the container height is chosen to he H — h while 
estimating the collision time with the lid. After colliding 
with the container, the particle may start to detach im¬ 
mediately from or stick to the container, depending on 
whether the criterion a > ai for the lid or a < Ub for 
the bottom is fulfilled or not. If not, the particle moves 
together with the container until the criterion is satis¬ 
fied. Subsequently the rupture process starts and the 
particle continues to bounce in the container. Based on 
the number of vibration periods the particle takes until 
it sticks back to the container bottom, we estimate the 
preferred periodicity for the wet granular layer under a 
certain vibration frequency and peak acceleration. 

As shown in Fig.[6l the numerical analysis predicts two 
classes of periodicity regions: One stays constant with / 
and the periodicity increases step by step with T, which 
we call periodicity class 1. The other one, called period¬ 
icity class II, starts at a certain frequency fc and grows 
with / monotonically, and the periodicity decreases as 
r grows. The spaces between neighboring periodicity 
regions are chaotic. This is manifested by the large ab¬ 
solute values as well as fluctuations of the periodicity 
arising from the immediate detachment of the particle 
after colliding with the container. 


Qualitatively, the class II periodicity regions agree with 
the pattern forming regions in the stability diagrams 
shown in Fig. [2] and 0] on the following features: (i) There 
exists a certain fc above which a certain periodicity re¬ 
gion starts, (ii) Both boundaries of a certain periodic¬ 
ity region grow with /. (iii) The period tripling region 
dominates and expands as / grows, (iv) The period dou¬ 
bling region starts at a lower fc and grows in a T region 
larger than the period tripling one. Different from the ex¬ 
periments, the model predicts an additional region with 
periodicity four below the period tripling one. 

In order to have a more quantitative understanding 
of the periodicity diagram and its dependency on the 
control parameters, we have varied both m and H the 
same way as the experiments. As a starting point, we 
explain why the periodicity is independent of / in peri¬ 
odicity class I. Following the above description, the gran¬ 
ular layer starts to detach from the container bottom as 
the acceleration of the plate satisfies a < a\) = —{n^ l)g 
in the laboratory frame. Assuming that the free flight 
period starts immediately (i.e., neglecting the rupture 
process), we can estimate the velocity at departure with 



where zq corresponds to the peak vibration amplitude. 
Using the definition of F and Ub, we have 


«c = ^yr2-(n + l)2, (3) 

which shows its dependency on the control parameters /, 
r and n. Note that the dependence on n represents the 
influence from sample mass m. As a first approximation, 
we further ignore the influence of vibration amplitude on 
the free flying period of the particle [i.e., assuming zq 
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FIG. 7. (Color online) Boundaries for period doubling (a) and tripling (b) regions from the numerical analysis. Various types 
of connection lines correspond to various combinations of parameters m and H, which are the same as in the experiments. 
Comparisons between the predictions of the numerical analysis [connected symbols, same legend as in panels (a) and (b)] and 
the onsets for both period doubling (c) and tripling (d) patterns measured experimentally in the F-F plane. See the text for 
a definition of the dimensionless control parameter F. 


{H — h)]. Consequently, the periodicity of the parabolic 
flight T can be estimated with 

07; 1 ,_ 

r=r— ■/l = r-x/r2-(n + i)2i, (4) 

g n 

which is independent of /, and grows with F as well as 
n. If we take the parameters used in Fig. [6] and skip the 
ceiling function [], Eq.|4] leads to F ^ 9.1 and 11.5 for 
periodicity T = 2 and 3, which fall into the middle of the 
ranges predicted by the numerical analysis (8.2-10.2 for 
periodicity 2 and 10.8-12.8 for periodicity 3). Thus, the 
periodicity class I corresponds to the situation without 
collisions with the lid of the container. 

As the above analysis applies as well to the cohesion¬ 
less case with n = 0, Eq. |4] gives us a clue on why the 
onset of subharmonic patterns observed in dry granu¬ 
lar matter [42| is weakly dependent on the vibration fre¬ 
quency. However, no instabilities are observed in this 
class of periodicity for the cohesive case that we focus on 
here, presumably due to the strong cohesion between ad¬ 
jacent particles as well as between the granular layer and 


the container. The granular layer may follow the period¬ 
icity collectively, but the energy injection is not sufficient 
to effectively mobilize the particles and to initiate insta¬ 
bilities. 

Eollowing a similar approach, we can understand how 
the lid plays a role in determining the periodicity class 
II. Taking the same assumption of Zo ^ {H - h), 
the criteria for the granular layer to reach the lid is 
Vc > y/2g{H — h). Together with Eq. [H we obtain the 
corresponding 


F = 


l87r^f^{H-h) 


-h (n + 1)2, 


(5) 


at which the particle starts to reach the lid. Erom there 
on, we expect the periodicity to decrease with the growth 
of T, because the time for the particle to travel within 
the container reduces. Qualitatively, Eq. [5] explains the 
growth of the threshold acceleration F with / and H ob¬ 
served in the experiments. Quantitatively, we cannot es¬ 
timate the F-f relation for various periodicities the same 
way as in the other case, because the periodicity relies not 
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only on the free flying periods but also on the time that 
the particle travels together with the container. 

Focusing on periodicity class II, we compare the pe¬ 
riodicity regions predicted for various m and H by the 
numerical analysis in Fig[71^a). It shows that either in¬ 
creasing H or decreasing m lead to enhanced onsets for 
both period doubling and tripling regions. This is in 
agreement with Eq. [5] (note that decreasing m is equiva¬ 
lent to increasing n). 

Moreover, Eq. [5] suggests a new parameter F = — 

h)/g that combines the influence of both / and H — h. 
It corresponds to a dimensionless acceleration where the 
amplitude is replaced with the maximum free travel¬ 
ing distance H — h of the granular layer. As Eig. [Tl^b) 
demonstrates, a replot of the numerical results shown in 
Eig. [Tfa) in the F-F plane leads to a collapse of the data 
points for a certain periodicity. Similar data collapse 
is found for the onset of patterns measured experimen¬ 
tally, except for the case of m = 81.6 g at relatively large 
F. The data collapse demonstrates the coupling between 
various control parameters (m, / and H) concerning the 
onsets of various periodicities in the system (note the 
relation between h and m shown above). 

A comparison between the numerical analysis and the 
experiments shows that the onset of period tripling pat¬ 
terns takes place at the upper bound of the periodicity 
region predicted by the numerics. This suggests that the 
periodicity is only a precondition for the patterns to ap¬ 
pear. Eollowing the above argument on why there are no 
patterns observed in periodicity class I regime, we specu¬ 
late that sufficient energy injection is necessary to gener¬ 
ate kinks, which are the basic ingredient for the standing 
waves as well as rotating spiral patterns observed exper¬ 
imentally. 

Moreover, good agreements between the experimental 
and numerical results on the critical F^ above which a cer¬ 
tain periodicity region emerges are found in Eig.[71^b). F^ 
can be estimated quantitatively as follows. Erom Eq. 01 
we estimate the acceleration to reach a certain periodicity 
T at a fixed n 


r « y7r2T2 + (n + l) 2 , (6) 

Supposing this is the acceleration at which the granular 
layer starts to touch the lid and the periodicity starts 
to decrease, we obtain, together with Eq. O the critical 
frequency fc at which the transition to periodicity T — 1 
starts: 


At the scale of the control parameter F, we estimate Fc 
for periodicity 2 or 3 to be T^/8 = 1.1 or 2.0, which 


agrees fairly well with the numerical (1.0 or 1.8) and the 
experimental (1.1 or 1.9) results shown in Eig. [7|^b). 

V. CONCLUSIONS 

In conclusion, the collective behavior of a thin layer of 
cohesive granular matter under vertical vibrations is in¬ 
vestigated experimentally. Stability diagrams for stand¬ 
ing wave patterns, rotating spirals and phase separations 
(granular ‘gas bubbles’) are presented. Eor the period 
doubling induced standing wave patterns, the granular 
layer splits into two domains that collide with the con¬ 
tainer at different phases. The interface between the do¬ 
mains corresponds to kink fronts that swing back and 
forth with time, in contrast to the case of period tripling 
induced rotating spirals . A comparison between the 
stability diagrams measured with various parameters re¬ 
veals that the phase separation region expands as the 
sample mass m (or correspondingly granular layer thick¬ 
ness h) decreases or the container height H increases. 
Moreover, the influence of these parameters on various 
pattern forming regions is discussed through comparisons 
between various stability diagrams. 

Using a simplified model considering the whole granu¬ 
lar layer as a single particle colliding completely inelas- 
tically with the container, we predict various periodicity 
regions and compare with the stability diagram measured 
experimentally. Eor both standing wave and rotating 
spiral patterns, we find corresponding period doubling 
and tripling regions in the periodicity diagram predicted 
by the model. Eurther comparisons indicate that addi¬ 
tional energy injection through colliding with the lid of 
the container is a precondition for both period doubling 
and tripling patterns to appear. Quantitatively, we in¬ 
troduce a dimensionless parameter F that characterizes 
the dependency of the boundaries for various periodic¬ 
ity regions on the vibration frequency /, height of the 
granular layer h as well as height of the container H. 
Moreover, the critical frequency fc above which a certain 
periodicity region starts can also be rationalized with the 
model. 

According to the model, the system could in principle 
generate patterns with a controlled number of periodic¬ 
ity (e.g., four-phase patterns jisMisj f that may lead to 
novel pattern forming scenarios. Eollowing the predic¬ 
tions from the model, a variation of the control parame¬ 
ters to explore possible period four instabilities will be a 
focus of further investigations. Moreover, we would also 
like to explore the possibility to control the dynamics of 
the domain interfaces with biharmonic driving miii. 
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